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Abstract 

We calculate the elliptic flow parameter V2 for Au+Au collisions at y/s~NN = 200 GeV employing 
the parton cascade BAMPS (Boltzmann Approach of Multiparton Scatterings). Besides gluon 
interactions including the bremsstrahlung process, interactions with quarks are considered in an 
effective, but approximate way to investigate the dependence of the collective flow on the number 
of parton degrees of freedom. We find that vi as a function of the transverse momentum pj- is 
sensitive to the number of parton degrees of freedom, whereas the pt averaged V2 does not. When 
including quarks, V2(pt) shifts to lower pr, the parton transverse momentum spectra become softer 
and the mean parton transverse momenta decrease. 
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I. INTRODUCTION 



Measurements on collective flow in heavy ion collisions at the BNL Relativistic Heavy 
Ion Collider (RHIC) [l| showed that the matter produced, the quark gluon plasma (QGP), 
is a nearly perfect fluid 2], Because the strength of the collective flow increases with 
decreasing viscosity, indirect extractions of the shear viscosity have recently been performed 
w tuning the shear viscosityfor the QCD coupling) as a parameter in viscous hydrodynamic 
4H8| and transport mo dels ;o match the elliptic flow V2- A consistent result has been 

achieved: on a conservative basis the shear viscosity to the entropy density ratio r]/s of 
the QCD matter at RHIC is less than 0.4 12]. Uncertainties stem from assumptions in 
modelling the various stages that the matter undergoes during its evolution. 

The results from the transport calculations using the Boltzmann Approach of Multiparton 
Scatterings (BAMPS) UJ showed that ^(pt) as a function of the transverse momentum p? 
is lower and the spectra are harder than the experimental data, whereas the integral of 
both gives a pt averaged t>2, which matches the experimental data. This inconsistency may 
stem from the nature of the rather complicated hadronization process than the simple parton- 
hadron duality picture used in llj. Another reason may lie in the fact that only gluons 
are considered as interacting constituents in the calculations presented in 11]. Including 
quark dynamics will increase the number of parton degrees of freedom and, thus, may soften 
the pt spectra and enhance t>2(j>r)- I n this paper we employ BAMPS and include effective 
quark degrees of freedom to investigate how the increase and the equilibration of partonic 
multiplicities affect the elliptic flow and also the 7]/s ratio of the QGP. 

This study shall demonstrate whether the buildup of elliptic flow in ultrarelativistic heavy 
ion collisions depends on the chemical equilibration of gluons and quarks. The latter de- 
termines the actual number of constituents in the partonic phase. In hydrodynamic calcu- 
lations it is assumed that the matter stays in chemical equilibrium and thus the number 
of constituents is conserved. However, the situation in a real heavy ion collision might be 



considerably more complicated 13 



14 



rium (if this occurs) later than gluons 



. Quarks are expected to achieve chemical equilib- 
131 ]. Even for gluons only, the chemical equilibration 



will proceed faster at the collision center than at the region near transverse edge of the 
parton system. The parton chemical equilibration can well be studied using the microscopic 
transport model with multiple scatterings such as BAMPS [15| . because BAMPS imple- 
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merits perturbative QCD (pQCD) elastic collisions as well as the pQCD based inelastic 
bremsstrahlung incorporating full detailed balance. 

The paper is organized as follows. In Sec. [Ill we review the parton cascade BAMPS and 
introduce effective quarks into BAMPS as a simplification for real quark dynamics. In Sec. 
IIHI the numerical results are shown to demonstrate the effect of the parton multiplicities on 
the elliptic flow of the parton matter. We summarize in Sec. IIVI 



II. BAMPS INCLUDING EFFECTIVE QUARK DEGREES OF FREEDOM 



The detailed model description of the on-shell parton cascade BAMPS can be found in 
Refs. 3,Q|. In short, the feature of BAMPS is the successful implementation of particle 
number changing processes with full detailed balance. This is ensured by using the stochastic 



interpretation of the transition rates in t 
cascade approaches can be found in [10, 



re Bo 
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tzmann equations for partons. Other parton 



The cross section of pQCD gluon elastic scatterings is given by 

da 99 ^ 99 9na 2 
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and the effective matrix e 



Gunion-Bertsch form 13, 



ement of pQCD inspired bremsstrahlung gg -h- ggg is taken in a 
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(3) 

where g 2 = 47ra s and do = 16 is the gluon degeneracy factor for N c = 3. denotes 
the perpendicular component of the momentum transfer, (k± = |kj_|) the perpendicular 
component of the radiated gluon momentum and y its rapidity in the center-of-mass frame 
of the collision, respectively. The suppression of the bremsstrahlung due to the Landau- 
Pomeranchuk-Migdal effect is effectively taken into account within the Bethe-Heitler regime 
using the step function in Eq. ([2]). Gluon radiations and absorptions are only allowed if the 
formation time of the process, typically r = coshy/k±, is shorter than the mean free path 
A m fp of the radiated or absorbed gluon. A m f p is calculated self-consistently [15|, |24j, i.m., 
A m fp is the inverse of the total collision rate R gg -^ gg + R gg ^ ggg + R ggg ^ gg , where R gg ^ggg and 
R g g g ^gg also depend on A m f p as indicated in Eqs. fl2]) and ()3]) [see Eqs. ([8]), ([9]), ( Tl^j) . and 

(USD]- 



The pQCD interactions are Debye screened. The screening mass is given by [1 31 ] 

/dp 1 
j—^-[N c f g (x,t,p)+n f f q (x,t,p)], (4) 

where n/ is the number of quark flavor, and f g and f q are the distribution function of gluons 
and quarks at a certain quantum state. For a pure gluon matter f q = 0. 

The initial parton distribution for BAMPS calculations is the same as that chosen in 
Refs. 0, [llj|: an ensemble of gluon minijets with transverse momenta greater than po = 1.4 
GeV, produced via semihard nucleon-nucleon collisions within a Glauber picture. The value 
of po was chosen by matching the parton cascade result of the final transverse energy per 



rapidity to the experimental data [9 



11] . Quark minijets, which take about 20% of the total 



parton number, are neglected for present studies, in order to easily make comparisons with 
the previous results obtained in [ll|. 

To take into account the quark-gluon changing process gg -H- qq and the elastic as well as 
the bremsstrahlung process involving light quarks, we effectively enlarge the gluon degener- 
acy factor d G . With a larger d G than 16, the ggg — > gg process is suppressed [see Eq. (J3J)] 
and thus more particles (i.e. quarks) will be produced. This numerical prescription indicates 
two assumptions: (1) Quarks are massless. (2) The elastic and the bremsstrahlung process 
involving quarks are identical with those of gluons according to Eqs. (JT|)- (J3j) . These assump- 
tions lead to the same kinetic equilibration of quarks and gluons, i.e., f q ~ f g . The second 
assumption overestimates the interaction rate of a quark approximately by the color-charge 
factor 9/4 due to the dominance of gg —> gg, qg — » qg, gg -H- ggg and qg <-> qgg in kinetic 
equilibration. The real equilibration of quarks is expected to be slower than that of gluons. 
Moreover, the present prescription used to include quarks cannot separate quarks from glu- 
ons. The relative fraction of the quark (or gluon) number to the total parton number is not 
known. Quark thermalization by real quark dynamics will be investigated in the future 25 ]. 
The simplified version used here serves as a tool to study the effect on the elliptic flow when 
varying the number of parton degrees of freedom. 

The enlarged do value can be obtained when considering a fully thermalized system of 
quarks and gluons. In this case f q = f g and the parton distribution function is given by 
/ = dcf g + dqf q = (do + dg)f g , where g?q = 24 is the quark degeneracy factor for n/ = 2. 
To include effective quark degrees of freedom in BAMPS we enlarge dc from 16 for a pure 
gluon system to 40 for a quark gluon system with two flavors. 
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Assuming f q = f g , which is also valid when the kinetic and chemical equilibration of 
quarks and gluons proceeds identically, the screening mass Eq. (j4]) becomes 

2 r d 3 p 1 

™d( x ^) = 3 nasN cj 7^33 - /( X '*'P) > ( 5 ) 

which is a factor of 2/3 smaller than the value at the beginning of the expansion, because 
initially there are only gluons. To make a reasonable description of the early stage, we use 
instead of Eq. (JSJ) 

/d 3 p 1 
J^y ~ /(X,t,p) (6) 

during the entire parton evolution. Accordingly, m 2 D is overestimated by a factor of 1.5 at 
the late stage of the expansion, when partons thermalize. The true screening mass will be 
changing from Eq. ([6]) to Eq. (jSJ) in time according to the true chemical equilibration of 
gluons and quarks, which, however, cannot be demonstrated in the present studies. More 
discussions will be given in the next section. 

As already considered in Refs. for the present BAMPS calculations, the kinetic 

freeze-out of particles occurs when the local energy density drops below e c , which is assumed 
to be the critical value for the occurrence of hadronization. In this paper we set e c = 1 
GeV fm~ 3 , which leads to a critical temperature T c = 200 MeV for a p ure gluon plasma 
and T c = 160 MeV for a quark gluon plasma with two quark flavors 26J. After the freeze- 
out partons are regarded as massless pions according to a simple parton-hadron duality 
picture. To consider realistic chemical and kinetic freeze-out a hadronization model and the 
subsequent hadron cascade should be included to BAMPS, which will be done in the future. 
This, of course, can also have certain influence on the findings in the next section. 



III. RESULTS 



With the assumptions for the quark dynamics in BAMPS we calculate the space time 
evolution of the quark and gluon matter produced in Au+Au collisions at tJsnn = 200 
GeV. The coupling is set to be a constant a s = 0.6. We evaluate the elliptic flow parameter 
t>2 as the average of (p x — Py)/{p x +Py) over particles within a certain window of momentum 
rapidity y = | ha[(E + p z )/(E — p z )]. Figure [T] shows the buildup of the elliptic flow v 2 at 
midrapidity \y\ < 1 in a Au+Au collision with an impact parameter of b = 8.6 fm. The 
solid curve gives the result of a pure gluon matter, which was already obtained from our 
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FIG. 1: (Color online) Time evolution of the px averaged t>2 for particles within |y| < 1. 

previous work [9]. The dashed curve shows the new result including quarks effectively with 
the enlarged dc = 40. We see a perfect agreement between the two results except for a 
difference of about 5% at later times. The key quantity to understand this agreement was 
oroposed to be the collision rate per particle R = n{a) (or similarly the Knudsen number 
271]): For a fixed collision geometry and a given initial condition the t>2 generation depends 
only on R, but not on the particular details of interactions among constituent particles. 
With more degrees of freedom (larger do) the total particle density n becomes larger due to 
the production of quarks. On the other hand, the screening mass ([6]) also becomes larger, 
which decreases the total cross section {a) of all the interaction types [see Eqs. ([I])-©]. The 
increase of n and the decrease of (cr) could keep the collision rate per particle and thus the 
t> 2 unchanged. 

To obtain a quantitative answer to this issue, we calculate numerically the total collision 



rate per particle R gg ^ gg + R gg ^ ggg + Rggg-> gg , where 



R, 



R, 
R 



99^99 



99^999 



n(v re \a gg ^ gg )2 
n(v re \a gg ^ ggg ) 2 



999^99 



-n 



999 ^99 



8E1E2E; 



3/ 3 



(7) 
(8) 
(9) 
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The averages are defined as 

(0) 2 = 1 J dT^Ti h(x, Pl )f 2 (x,p 2 ) O, (10) 

(Q) 3 = dV x dV 2 dV 3 fi(x, Pl )f 2 {x,p 2 )f 3 {x,p 3 ) Q (11) 

with 

n = JdT 1 f(x,pi), (12) 

where dTj = d 3 pj/(2ir) 3 , j = 1,2,3, and the local particle distribution function f(x,p) is 
the output from the parton cascade, f(x,p) = Y,i$ [x — Xi(t)} 5^[p — pi}. The sum runs 
over all particles with individual position Xi and momentum Pi at time t. 

For a gg — » gg or a gg — >■ ggg process involving two incoming particles with (x,pi) 
and (x,p 2 ) the relative velocity of the two particles is given by v Ie \ = s/(2EiE 2 ), where 
s = (pi + P2) 2 is the invariant mass. The total cross sections read 

1 r s / 4 . n da„ 



1 j « C /r;|^ M9 | 2 (2vr)V 4 )(p 1 +p 2 - pi -ft - ft) , (14) 
where gT^- = d 3 p'j / (2tt) 3 / (2Ej) , j = 1, 2, 3. For a ggg — )■ ^ process we define 

1 999->99 =2\j dT 'l dT '2\ M 999^99? ( 27X Y 6{4 \Pl + P2 + P3 ~ Pi ~ P 2 ) , (15) 

which is, similar as a cross section, an integral over the final states. The factor of 1/2! in Eqs. 
(TIB"]) and (ITS]) is due to the fact that the two outgoing gluons are identical particles. There 
is no such a factor (1/3!) in gg — > ggg processes [28(, because the three outgoing gluons 
are distinguished as the radiated gluon (p 3 ) : the radiating gluon (p 2 ) and the gluon ( Pl ), 
which is only deflected and does not radiate. Thus, they are kinematically "distinguishable" 
particles. 

We note that n in Eq. (fT2|) is the local particle density in the center of mass frame of the 
Au+Au collision (collision frame) and only equals the Lorentz invariant particle density if 
the local particle system is at rest, e.g. at the collision center. Also the local collision rates 
(JIJ)-(J2J) are calculated in the collision frame and are smaller than those in the rest frame by 
a Lorentz factor. 

The left panel of Fig. [2] shows the time evolution of the total collision rate per particle 
{Rgg^gg + Rgg^ggg + Rggg^gg) . The rate is obtained as the average over all particles in 
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FIG. 2: (Color online) Time evolution of the collision rate (left) and the transport collision rate 
(right) per particle at midrapidity \rj s \ < 0.2. 

the whole transverse plan at the central space time rapidity \r] s \ < 0.2. r] s is defined as 
T] s = | ln[(t + z)/(t — z)\. The local rate per particle is calculated according to (JTJ)- (|9j) . We 
have also computed the average rate by counting interaction events, which occur in BAMPS 
within certain time intervals. The results from both calculations agree with each other. 

As mentioned before, the collision rate per particle, may determine the buildup of the 
elliptic flow t> 2 , which is shown in Fig. [1] within the momentum rapidity \y\ < 1. We have 
also calculated t>2 within \q s \ < 0.2. The agreement of the two curves seen in Fig. [T] is 
unchanged. The reason why we show the v-i result within \y\ < 1 is to compare it with the 
experimental data (see Fig. [3]). 

From Fig. [2] we see a moderate difference between the collision rates in the simulations 
with do = 16 and do = 40, whereas the average transport collision rates per particle shown 
in the right panel are much closer to each other. The slightly smaller transport collision 
rate in the calculation with dc = 40, compared to that with do = 16, is the reason for the 
slightly smaller final t>2 with do = 40 than that with do = 16 (see Fig. [1]). 
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The local transport collision rate per particle is given by 16] 



R 



tr 



fdT^Qlfj-i^Jd^Qlf} 



71 ' 3 



(16) 



where (O) = (l/n) JdT 1 fiix^o^O and d[f], i = gg ->■ g^,^ ->■ ggg,ggg ->■ denotes 
the respective collision term 16 L The transport collision rate per particle quantifies the time 
scale of thermal equilibration 16]. Because the angle #1 in p 2 zl /E\ = cos 2 #1 enters both the 
gain and loss term of Cj[/], the larger the mean momentum deflection in each collision, the 
larger is the transport collision rate and the faster is the momentum isotropization. The 
latter leads to a faster buildup of the pressure gradient and thus to a larger elliptic flow. 
Therefore, the transport collision rate, rather than the collision rate, relates more closely to 
the buildup of the elliptic flow V2, as demonstrated in Fig. [2] in relation with Fig. [U 

The resulting smaller difference in the total transport collision rates than that in the total 
collision rates (when varying da from 16 to 40) stems from the time evolution of the screening 
mass m D [see Eq. 



which determines not only the cross sections a gg ^g g , cr g g-^gg g , and 



1 ggg^gg 



, but also the collision angle in these processes. As mentioned after Eq. 



the 



choice of mo is physically reasonable, because the dynamics at the early stage is dominated 
by gluons. On the other hand, is overestimated at the late times, when "quarks" have 
considerable fraction of the whole parton system. Thus, the cross sections and the rates at 
the late times shown in Fig. [2] are somewhat smaller than their true values, when the fully 
realistic quark dynamics is considered. 

In Fig. [3] we show the final elliptic flow t>2 within \y\ < 1 for various mean numbers of 
participating nucleons (N pai t) (equivalently impact parameters b [ill). The open squares 
(or the solid curves) in Fig. |3] and all the following figures depict the results from the 
calculations with do = 16, which have already been presented in Ref. llj. The open circles 
show the new results with da = 40, which are slightly smaller than the values calculated 
with do = 16 (as also seen in Fig. [T|), and are still comparable with the experimental data 



at RHIC 
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Figure H] presents the final transverse energies at midrapidity. If quarks are included in 
the initial condition for the parton cascade calculations, the results would be closer to the 

n 

experimental data [31]. Comparing the results obtained in the calculations with da = 40 
to those with da = 16, we find a tiny difference. This implies the same decrease of the 
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FIG. 3: (Color online) Elliptic flow parameter v% (pr averaged) at midrapidity as a function of 
number of participating nucleons. 
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FIG. 4: (Color online) Transverse energy at midrapidity as a function of number of participating 
nucleons. 
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FIG. 5: (Color online) Mean transverse momentum at midrapidity as a function of number of 



participating nucleons, compared with the STAR data 
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transverse energy during the expansion, which is the consequence of the same viscous effect 
regardless of the details of microscopical processes. 

On the contrary, the mean parton transverse momentum shown in Fig. [5] does depend on 
the details of microscopical processes. The larger the number of parton degrees of freedom, 
the smaller is the final mean transverse momentum. Suppose the system is in full thermal 
equilibrium, we obtain the energy density e ~ dcT 4 and the particle density n ~ dcT 3 ~ 
dc^e 3 / 4 . Similar as the decrease of the transverse energy the local energy density is not 
sensitive to da- Thus, n ~ d 1 ^ and the particle number with do = 40 is a factor of 1.26 
more than that with da = 16. Therefore, to the maximum (in case of thermal equilibrium) 
the mean transverse momentum will be reduced by a factor of 1.26, when do is enlarged from 
16 to 40. From Fig. [5] we realize an average reduction by a factor of 1.2, which is less than 
the maximum reduction. This indicates that full thermalization is not immediately achieved 
and particularly the net particle production (chemical equilibration) in the calculations with 
do = 40 is less complete at the freeze-out than for the pure gluon system (<i G = 16). The 
reason is that for the same initial condition the system with do = 40 is initially farther 
apart from full chemical equilibrium than the system with da = 16 (see Fig. [TOi) . Because 
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FIG. 6: (Color online) Momentum spectrum for central collisions. 

the production rates R gg ^gg g are almost the same in both cases, the chemical equilibration 
of the system with da = 40 proceeds always behind the process in the system with dc = 16. 

The production of more particles from the same energy content leads to a smaller tem- 
perature at the (chemical) freeze-out. This clearly relates the softening of the transverse 
spectrum dN/ {pxdpTdy) shown in Fig. [6] as an example for most central collisions. Both re- 
sults of the mean transverse momenta and the transverse momentum spectrum come closer 



to the experimental data 3jJ, l32| , if the effective quark degrees of freedom are included into 
the parton cascade calculations. 

On the one hand, the pt averaged elliptic flow v% is almost unchanged, when enlarging da 
from 16 to 40. On the other hand, the slope of the transverse momentum spectrum increases 
due to the production of more particles with larger da- Thus, there must be a simultaneous 
change in the differential elliptic flow ^(pt), because t>2 = / dpr ^(pt) dN/N/dpr- The 
change of i>2(pr) is shown in Fig. [7] for the most central 50% collisions and in Fig. [8] for the 
mid-central class. The dotted curves are just the dashed curves times an effective factor, 
which enhances the pt averaged v 2 values of the calculations with do = 40 to be equal with 
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FIG. 7: (Color online) Momentum dependence of the elliptic flow V2(pt) for the most central 50% 
collisions. 
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FIG. 8: (Color online) Momentum dependence of the elliptic flow V2(pt) for the centrality bin 
30 - 40%. 
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those calculated using do = 16. The comparison between the dotted and the solid curve is 
then equivalent to that of the i^Cpt)/^ scaling. Unlike the px averaged t>2, ^(pt) depends 
on the number of parton degrees of freedom: on average, ^(Pt) moves up 20% or shifts 20% 
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to the left toward small momentum and becomes closer to the experimental data 
when quarks are added. This is the major finding of this work. 

V2(pt) values at large pt decreases by 20% and is a factor of two smaller than the experi- 
mental data. Although this has a negligible effect on the pt averaged t>2, because the parton 
number there is tiny, as seen in Fig. [6j it is important to understand hadronization. From 
our results we see a bend at pr ~ 1-2 GeV. For higher t>2(i>r) saturates. If parton re- 



combination models [34J govern the hadronization of partons with high p?, i^Qot) of mesons 
will be twice of partons. This gives a reasonable explanation of the factor of two difference 
between the experimental data and the parton cascade v-i results at high pt- On the other 



hand, viscous hydrodynamic calculations [6] showed that v 2 at high p T is very sensitive to 
the value of shear viscosity. Also, the difference of the viscous correction to the equilibrium 
phase space distribution between mesons and baryons at freeze-out can explain constituent 
quark scaling without quark recombinations [35]. Hadronization is still an open issue. 
As a final result the shear viscosity rj from the parton cascade simulations is extracted 



by using the formula, which has been derived in 36] and applied in |9j: 



i , 17 , 

71 ~ <=, i — / 7 ,2\ y\ E>tr I 1 k c> _d ' \ L( ) 

3 \ u zl l^i ±L i ~ L - Ull 'gg^ggg 1L ggg^-gg 

where i = gg — >■ gg,gg — > ggg,ggg — > gg arid v z = p z /E. The transport collision rate 

R\ T is obtained via Eq. (JTBjl . The entropy density s is calculated by assuming the kinetic 

equilibrium, which gives 

s = 4n-nlnA, (18) 

where A = n/n eq with n eq = dcT^/it 2 and T = e/(3n) defines the gluon fugacity. The 
true entropy density is expected to be slightly smaller, because overall kinetic equilibration 
cannot be complete in an expanding viscous system. 

Figure IH] shows the rj/s ratio in the center of a Au+Au collision at the impact parameter 
b = 8.6 fm. Results from two calculations with do = 16 (solid curve) and da = 40 (dashed 
curve) are compared. At early times the rj/s ratio from the calculation with da = 40 is 
smaller than that with do = 16. Because the particle density n in rj and s cancels, and 
quantities R\ r as shown in Fig. [2] and (E(l/3 — i>f))/(l/3 — {v 2 z )) are not sensitive to the 
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FIG. 9: (Color online) Time evolution of the ratio of the shear viscosity to the entropy density. 

number of parton degrees of freedom, therefore, the difference in 7]/s comes mainly from the 
different chemical equilibration in two calculations, i.e., the difference in fugacity A, which 
appears in Eq. (fTT|) as 1.5R gg ^ ggg — R gg g^gg ~ 1-5(1 — X)R gg ^ ggg and in Eq. ffT8]) as In A. 

In Fig. [10] we see that the parton fugacity A(t) of the system with do = 40 is smaller 
than that with do = 16 due to the increase of parton degrees of freedom after the initial 
production. Smaller A decreases r\ and increases s, and, thus, decreases the rj/s ratio, 
especially at early times when the system is far from the chemical equilibrium. We note 
that the total entropy density of the quark and gluon system can only be calculated via 
Eq. (ITS]) , if the A values of quarks and gluons are the same. In the reality, gluons dominate 
the early stage of the plasma and the gluon fugacity should be larger than the quark's. 
Therefore, the true entropy density is smaller than that obtained via Eq. (flHl) . The true r]/s 
ratio of a quark gluon plasma should be larger than that presented as the dashed curve and 
thus may become closer to the values for a pure gluon plasma. Moreover, the contribution of 
the chemical equilibration to rj, the term 1.5R„ g ^ ggg — R ggg ^ gg in Eq. ( !T7|) . does not appear 
explicitly in the expression derived in Ref. |37| by using the second order Grad's expansion. 
This may slightly increase rj compared to the Navier-Stokes value via Eq. ( flTl) . if the system 
is out of chemical equilibrium. 
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FIG. 10: (Color online) Time evolution of fugacity. 

IV. SUMMARY 

Employing the on shell parton cascade BAMPS we have studied the effect of an increasing 
number of parton degrees of freedom on the elliptic flow parameter v<i generated in Au+Au 
collisions at the RHIC energy ^snn = 200 GeV. The initial condition for BAMPS is assumed 
to be an ensemble of gluon minijets. The additional effective quark degrees of freedom during 
the further evolution are created by increasing the degeneracy factor da from 16 for a pure 
gluon system to 40 for a quark gluon system. This prescription indicates the assumption that 
quarks and gluons are identical particles. With this assumption for the BAMPS calculations 
we have found that the pt averaged v 2 values and the total transverse energy at midrapidity 
are almost unchanged with or without quarks, which is a consequence of the almost same 
transport collision rates during the entire expansion. Second, incorporating quarks the 
parton multiplicities at freeze-out increase, which leads to a decrease of the mean parton 
transverse momenta (px) and a softening of the transverse spectra. Simultaneously, the 
differential elliptic flow ^(pr) shifts toward lower momentum. Adding quarks with two 
flavors brings a 20% effect on (p?) and V2(pt), which is smaller than the maximum value 
of 26% due to the incomplete chemical equilibration in the present studies. The incomplete 
chemical equilibration is also the reason for the slightly smaller rj/s value in the quark gluon 
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plasma compared with the result for a pure gluon plasma. 

The present prescription for the inclusion of quarks gives an estimate of interactions 
with quarks. Bremsstrahlung involving quarks should and will be explicitly implemented in 



BAMPS in the future 



251 ]. because this process is essential for quantifying the thermalization 



and the elliptic flow of quarks. Moreover, details on the elliptic flow of both quarks and 
gluons at the phase transition will provide a more quantitative basis for understanding the 
possible effect on the elliptic flow due to the hadronization of the deconfmed matter. 
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